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I. INTRODUCTION 

The one-body matter density distribution plays an im- 
portant role in our understanding of nuclear structure. In 
nuclei it determines the intranuclear distribution which 
includes both protons and neutrons, unlike the charge 
density which determines the proton distribution. It re- 
quires hadronic probes and for this reason it is not as ac- 
curately determined as the charge distribution. However 
at new experimental facilities such as RIKEN in Japan 
there are plans to improve the measurements of the mat- 
ter density distribution using proton beams Simi- 
larly matter density distributions of individual hadrons 
are determined in hadron proton collisions . From the 
analysis of elastic differential cross sections for n^p and 
K^p, the radii of these mesons can be extracted within 
certain assumptions like the eikonal approximation and 
the exponentiation of the S matrix 3] . 

In this work we study hadron matter densities using 
gauge invariant correlators calculated in lattice QCD. 
This is an extension of our previous work Q on hadron 
charge density distributions. In addition to the mat- 
ter density distribution we present results for the pseu- 
doscalar density and compare to bag model predictions. 
Both quantities are evaluated in quenched QCD and with 
two flavours of dynamical quarks. To test our lattice pro- 
cedure we perform a more detailed analysis of the rele- 
vant correlation functions as compared to what was done 
in ref. iJ] extending it also to the matter and pseudoscalar 
densities. By comparing the matter and charge density 
distributions one can draw important phenomenological 
conclusions regarding hadron deformation. 



II. GAUGE INVARIANT CORRELATION 
FUNCTIONS 



We consider the equal-time correlators 




C^{v,t) = j d\' {H\'j^{v',t)j^{v' + v,t)\H) , (1) 

with the current jp(r, t) given by the normal order prod- 
uct : u{Y,t)Tu{r,t) :. For F = 70 and F = 1 we obtain 
the charge and matter density distributions respectively. 
The pseudoscalar density obtained with F = 75 is also 
evaluated since it provides a useful observable for testing 
the predictions of the bag model. 

The matrix elements for mesons and baryons consid- 
ered in this work are shown schematically in Fig. ^ In 
the case of the baryons two relative distances are involved 
and three current insertions are required. However we 
may consider integrating over one relative distance to ob- 
tain the one-particle density distribution shown schemat- 
ically in the lower part of Fig. ^ Baryon matrix elements 
with three current insertions, one on each of the three 
quark lines, were computed in the case of the charge den- 
sity distribution in ref. ^ and shown, after integration 
over one relative distance, to reproduce, within statistics, 
the one-particle density. In this work we will only con- 
sider one-particle distributions, and therefore only dia- 
grams with two current insertions such as shown in Fig.^ 
are evaluated. In the future, we plan to extend our work 
to three current insertions, at equal and unequal times. 

We also address here an important technical point. 
On our lattice of time-extent T, with the usual anti- 
periodic boundary conditions for the fermion fields, we 
insert the currents at Euclidean time separation T / 4 from 
the source and from the sink, as indicated in Fig.^ The 
current insertions must be separated far enough from the 
source and the sink to suppress excited hadronic states 
and non-zero momenta. The latter are not projected out, 
as often done in other studies, because zero-momentum 
projection requires a summation over spatial translations 
of the source or sink. Here, this summation is technically 
not feasible because it would involve quark propagators 
from all to all spatial lattice sites. In the following sec- 
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FIG. 1: Equal-time current-current correlator for a meson 
(upper) and a baryon (lower), t and T /2 — t must be large 
enough to filter out the hadronic state. T is the time exten- 
sion of our lattice. We have assumed anti-periodic boundary 
conditions in the temporal direction in drawing this figure. 



tion we will demonstrate that our time separation of T/4 
is sufficient to satisfactorily filter out undesired states. 
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FIG. 2: The local mass, m(t), at k = 0.153 (a) for the 
pion and (b) for the rho is shown versus the time separation 
from the source in lattice units, for a local source with anti- 
periodic boundary conditions (x), a local source with Dirich- 
let boundary conditions (crosses), a smeared source (stars) 
and a wall source (filled triangle). The dashed line shows 
the plateau value obtained by fitting the data coming from 
Dirichlet boundary conditions. 



III. LATTICE TECHNIQUES 



All the results presented in this work have been ob- 
tained on a lattice of size 16^ x 32. For the quenched 
case we use 220 NERSC 13 configurations generated at 
(3 = 6.0, and for the unquenched case we analyse two k 
values using for each 100 SESAM configurations ^ sim- 
ulated at (3 — 5.6 with two degenerate quark flavours. 
The physical volume of the lattice in the quenched and 
in the unquenched case is approximately the same. 

Let us first demonstrate that the time extent of our lat- 
tice is large enough to isolate the hadron ground state of 
zero momentum. For this study we analysed 56 quenched 
configurations at k, — 0.153. The first check is to eval- 
uate the mass of the particles via the two point corre- 
lator, C(t), using different sources. For this comparison 
we use: i) a local source with anti-periodic boundary 
conditions in the temporal direction, ii) a local source 
with Dirichlet boundary conditions in the temporal di- 
rection, iii) Wuppertal smearing of the source and iv) a 
wall source which projects out the zero momentum quark 
propagator. In the latter case we Coulomb gauge fix the 
configuration at the time slice of the source. Employ- 
ing Dirichlet boundary conditions allow us to utilize the 



whole time extent of our lattice instead of using half as is 
done for anti-periodic boundary conditions in the rest of 
the cases. When using Dirichlet boundary conditions the 
source is placed at the second time slice, the sink at the 
last but one and the current insertions at various inter- 
mediate time separations, t, from the source. In all cases 
we use a local sink. The results for the effective mass, 
m(t), for the pion and the rho are shown in Fig. |21 As 
expected Wuppertal smearing of the source produces the 
earliest plateau. Local sources reach the same plateau 
value for time separations < > 8a from the source, where 
a is the lattice spacing. Using a wall source does not im- 
prove projection to the ground state. In fact the effective 
mass still deviates at f = 8a and only converges to the 
plateau value at time separations t > 10a. For testing 
the mass plateaus using two point hadron correlators, we 
have summed over the sink spatial volume which projects 
to the zero momentum state in the standard way. Had we 
not performed the sum over the sink spatial volume the 
errors on the effective mass would have be much larger 
and the comparison of different sources would have been 
less meaningful. As we have already pointed out sum- 
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ming over the spatial volume of the source is not pos- 
sible for the density distributions since it would require 
the all-to-all propagator. Therefore it is crucial to check 
that higher momenta are sufficiently suppressed for the 
evaluation of the density distributions. 
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FIG. 4: The charge density distributions at k = 0.153 versus 
the radius in lattice units (a) for the nucleon and (b) for the 
using a local source and sink with anti-periodic boundary 
conditions (x's), for a smeared source and sink (stars) and for 
a wall source and sink (filled triangle). Data obtained using 
Dirichlet boundary conditions are too noisy and they are not 
shown. 



FIG. 3: The charge density distribution at k = 0.153 (a) for 
the pion and (b) for the rho is shown versus the radius in lat- 
tice units when the current couples to the quark at time sepa- 
ration from the source of t = 8a with anti-periodic boundary 
conditions and oi t = 14a for Dirichlet boundary conditions. 
In the case of the rho distribution the errors on the data 
obtained with Dirichlet boundary conditions are omitted for 
distances greater than r/a = 6 for clarity. The notation is as 
in Fig.H 



In Fig. O we show the results for the charge density 
distribution at different time slices for the pion and the 
rho, obtained using Dirichlet boundary conditions which 
allow us to probe larger time separations. Here, in con- 
trast to the two point correlator, C(f), the source and 
the sink are treated in the same way. For all the plots 
showing the charge density distribution as a function of 
r we have averaged the lattice results over bins of size 
0.28a. As it can be seen the results for the pion are the 
same at quarter lattice time separation t — 8a and at 
t = 14a. We also have agreement between the results 
obtained with local and smeared sources. For the rho 
the comparison is more difhcult since the results using 
Dirichlet boundary conditions are more noisy than for 



the pion. Within our statistics the charge density distri- 
butions obtained at time separation t = 8a and t — 14a 
are the same. Also, apart from small deviations at short 
distances, there is agreement of the results obtained us- 
ing Wuppertal smeared and local sources. This means 
that any contamination from high momentum states is 
small and, within our statistical uncertainties, it is per- 
missible to use anti-periodic boundary conditions with 
current insertions at t = 8a. This is fortunate, since it 
allows us to analyze standard full QCD configurations. 
It is also advantageous since the gauge noise is far less 
for time separations 8a. Using a wall source and sink 
produces a correlator which is in disagreement with the 
rest at time separation t = 8a. This means that there is 
still a sizable contamination from excited states, which 
is not surprising given that the effective mass for this 
source has not completely converged at t — 8a. In con- 
trast to the other sources, where the individual quark 
propagators can carry non-zero momentum, wall sources 
project on zero momentum for each quark propagator. 
If the time evolution is too short to filter the ground 
state then the quarks carry lower momentum than in the 
exact ground state. This would explain the fact that. 
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for a wall source, the charge density distribution is still 
broader than the rest ai t — 8a: the quarks carrying a 
smaller momentum lead to a slower decay. This effect is 
less visible in the pion, which indicates that in the ground 
state the quark relative momentum is smaller than in the 
rho. In Fig. 01 we compare the results obtained with the 
different sources for the nucleon and the A+ at time sep- 
aration t = 8a. Data with Dirichlet boundary conditions 
are not shown since they are too noisy. Again using a 
wall source and sink produces a correlator that decays 
more slowly than the rest, showing that projection to 
the ground state is less effective when the initial and final 
states are constructed with zero momentum quarks. On 
the other hand, there is agreement between results ob- 
tained with smeared and local sources apart from small 
deviations at short distances for the A''" channel. The in- 
dependence of the results from the different interpolating 
fields indicates that the ground state has been isolated 
sufhciently well. Therefore for the parameters used here, 
it is justified to use anti-periodic boundary conditions 
and local sources to study the density distributions. Us- 
ing Dirichlet boundary conditions to increase the time 
separation between the source and the sink produces re- 
sults that are more noisy and can not be distinguished 
from the ones obtained with anti-periodic boundary con- 
ditions. 



In this work we use lattices having similar lattice spac- 
ing a. In order to assess finite lattice spacing effects on 
these distributions we compare results using the contin- 
uum current and the lattice conserved current which dif- 
fer by 0{a) terms. The lattice conserved current is given 
by 

/ 

-i>f{x)il-jf,)U^'ix)^fix + fL)} (2) 

and it is symmetrized at lattice site x. Q/ is the charge of 
quark of flavor /. The results for the unnormalizcd den- 
sity correlator, shown in Fig. |31 demonstrate that finite- 
a effects lead to the wrong behaviour near the origin. 
The same conclusion was also reached in the study of the 
density distribution of charmonium where results using 
the Wilson Dirac operator were compared to those using 
clover and tadpole improved operators 0- In the rest 
of this work, unless otherwise stated, we will normalize 
the density correlators over the spatial volume and not 
over their value at the origin because of the potentially 
large finite a error at the origin. An improvement beyond 
the present work will be to use an O(a)-improved Dirac 
operator to eliminate finite-a effects at short distances. 
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FIG. 5: (a) The pion, (b) the rho, (c) the nucleon and (d) A+ 
charge density distribution versus the radius in lattice units 
using the local current (x's) and the lattice conserved current 
(crosses) . 



IV. MATTER DENSITY 

In this section we present results on the matter den- 
sity distribution using a local source and a local sink. 
As discussed in the previous section, our best compro- 
mise for the study of the density distribution is to fix the 
source and the sink for maximum separation ai ti = a 
and tf = 17a and insert uu and dd operators half-way be- 
tween dX t — ti = 8a. By examining the dependence of the 
matter correlator on the time separation t — ti at which 
the uu and dd operators are inserted, as done above for 
the charge density, we again conclude that a time interval 
oi \t — ti\ — \t ~ tf \ — 8a is sufficient to project to the 
ground state of the hadrons of interest. In this evaluation 
we use the same interpolating fields, described in detail in 
ref. 3, as for the density correlator, and the Wilson Dirac 
operator with hopping parameter k — 0.15, 0.153, 0.154 
and 0.155. The ratio of the pion mass to the rho mass 
at these values of n is 0.88, 0.84, 0.78 and 0.70 respec- 
tively. Using the standard definition of the naive quark 
mass, 2mqa — (1/k — 1/kc), where Kc — 0.1571 is the 
value of K at which the pion becomes massless, we find 
ruq - 300, 170, 130 and 85 MeV respectively. To obtain 
these values we used the string tension [l^ to set the 
physical scale, which gives for the inverse lattice spac- 
ing a-^ ~ 1.94 GeV(a = 0.103 fm). This choice makes 
direct contact to our previous study of charge density 
distributions Q. Another possibility is to use the rho 
(or the nucleon) mass in the chiral limit. Using the rho 
mass yields a-i=2.3 GeV (a = 0.087 fm) [ullll, with a 
systematic error coming from the choice of fitting range 



5 



and chiral extrapolation ansatz of about 10%, which is 
about twice as large as the statistical error. We note that 
Ref . observed stability of that scale with respect to an 
increase of the lattice size. In our discussion of quenched 
data we will use the value of a determined from the string 
tension. However, to compare the quenched with the un- 
quenched results we will use the value extracted from the 
rho mass in the chiral limit, since this determination is 
applicable both in the quenched and in the unquenched 
theory. In our discussion of the bag model results we will 
set the scale using the nucleon mass since, in that case, 
the nucleon mass is used to fix the bag model parameters. 
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FIG. 6: (a) The pion, (b) the rho, (c) the nucleon and (d) 
matter density distribution versus the radius for n = 0.15 (x), 
K, = 0.153 (crosses) and k — 0.154 (stars). For the pion we 
also include the results for k = 0.155 (filled triangles). The 
errors bars are not shown for clarity. 

In Fig. 1^ we consider the quark mass dependence of the 
matter distribution. The pion shows the strongest depen- 
dence on the quark mass, in contrast to the charge dis- 
tribution where the rho showed the largest variation 0. 
For quark masses in the range of 300-100 MeV investi- 
gated here baryon charge and matter distributions show 
essentially no variation with the quark mass. However, 
charge and matter distributions are quite different from 
each other. The comparison between both is shown in 
Fig. [Tlfor K = 0.153. We observe, between the pion and 
the rho but not between the nucleon and the A"*", a larger 
variation for the charge density as compared to the mat- 
ter density distribution. In fact the matter density dis- 
tribution is very similar for the four hadrons considered 
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FIG. 7; Comparison of the charge and matter density dis- 
tributions at K = 0.153. Left for the pion (filled triangles 
for the charge) and (stars for the matter) and the rho (x 
for the charge) and (open rhombus for the matter) distribu- 
tions. Right for the nucleon (filled triangles for the charge) 
and (stars for the matter) and the (x for the charge) and 
(open rhombus for the matter) distributions. 



here. In all cases the matter distribution decays faster 
than the charge density, an observation consistent with 
the results of Ref. 
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FIG. 8: Rho asymmetry at k = 0.154 (a) for the charge and 
(b) for the matter density distributions. The upper curve 
labeled by Z is Cp(0,0, z) and the lower curves labeled X 
and Y is C^{x, 0, 0) and Cp(0, y, 0) respectively. The z-axis is 
along the spin direction of the rho meson. 



In Fig. we show the correlators for the spin-0 state 
of the rho, for quark separations along the spin axis (z) 
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FIG. 9: Decomposition of the density-density correlator, C^^, 
and the matter density correlator, Cj, for the rho meson into 
angular momentum part L = and L = 2 at k = 0.153 (left) 
and K = 0.155 (right). 



and perpendicular to it (cc). We observe a z — x asym- 
metry in the case of the charge distribution Q whereas 
in the case of the matter distribution no statistically sig- 
nificant asymmetry is present. The rho deformation seen 
in Fig. |H| can be made more quantitative by analysing 
the correlators into a dominant i = and a suppressed 
L = 2 state El : 



<Pfc(0)|j;(r)it(0)|pfc(0) >=/o(r) 



•'la 



3r2 



/2M 
(3) 

where \pk{0) > is a zero momentum state with polar- 
ization k. The deformation, (5, determined from the 
quadrupolc moment, can be defined as 



, 3 (3z2 - r2 
= r^TT — 



(4) 



Writing for the symmetry axis R3 — ^ {z^) — Rt + e 
where — l/2(x^-|-y^) and assuming a small symmetry 
deviation e we obtain 



(5 — J— ^3 ~ Rt 
Rt R 



(5) 



where R^ = l/3(x^ +y'^ + z'^) is the mean square radius. 
Expanding the hadron wave function in L = and L = 2 



components as ip{r) = a(j)o{r)Po{cos6)+f3(t)2{r)P2{cos9), 
where Pl are Legendre Polynomials, we find for the de- 
formation 



3 ^ (0ok2|02) 



5 a {(t)o\r^\4>o) 



(6) 



where we have neglected terms. Taking 4'Q{r) ~ 
exp(— mpr) and 4'2{r) ~ r^exp(— TO2r), and again ne- 
glecting 0^ terms which appear when squaring the wave 
function to obtain the correlator so that /o(f ) — a^(/)§(r) 
and f2{j') = 2>a(3(j)Q(r)4>2{r), we fit the resulting expres- 
sion to the lattice data. As shown in Fig. for two k 
values, this ansatz provides a good description of the 
rho charge density distribution. From the fits we find 
at the four k values that we have analysed a non-zero ra- 
tio P/a ^ 0.02 with an error of about 60% which mainly 
arises from the poor determination of the cocfhcient of 
the L — 2 state. This leads to a deformation of the order 
of 1%. The errors in extracting 6, within this approach, 
are large since, in addition to the ratio P/a, this evalua- 
tion also involves mg and (mo + 7712) amplifying further 
the error on 6. A direct determination of the quadrupole 
moment via Eq. ^yields i5 = 0.03 ± 0.01, with a better 
control on the errors than the value obtained from the 
angular decomposition. Unquenching tends to increase 
this deformation, but the statistical error also increases, 
so that no definite conclusion regarding the importance 
of pion cloud contribution can be reached. However, the 
fact that deformation is seen for the rho meson in the 
quenched theory casts doubt on models which assume a 
spherical core and attribute the deformation entirely to 
the pion cloud . The same analysis for the matter den- 
sity distribution yields an L = 2 component consistent 
with zero, as expected from the lack oi z — x asymmetry 
demonstrated in Fig.|Sl Also in the unquenched case, the 
matter density asymmetry remains consistent with zero. 

The absence of a measurable deformation in the case 
of the rho matter distribution has important implications 
for the mechanism that produces the asymmetry in the 
rho charge distribution. Since the matter and charge op- 
erators have the same non-relativistic limit, this strongly 
suggests that hadron charge deformation is a relativistic 
effect. That is particularly interesting since it has strong 
implications for the validity of various models used in 
the study of nucleon deformation as well as in the eval- 
uation of transition matrix elements for — > A. In 
a number of these models 16J the quadrupole moments 
are evaluated in the non-relativistic constituent quark 
model limit with D-state admixture or with two-body 
currents where one would expect the same results for the 
charge and matter distributions. In these models one 
thus expects that for the rho meson the charge and mat- 
ter quadrupole moments are the same up to an overall 
multiplicative factor related to the different electromag- 
netic and strong coupling constants. The conclusion that 
the charge deformation is a relativistic effect is confirmed 
by another study employing heavy quarks on a fine lat- 
tice. For these heavy quark systems no signal for charge 
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FIG. 10; Comparison of the quenched matter density distri- 
bution at K = 0.153 (left) to the unquenched one at k = 0.156 
(right) for the pion (x's), the rho (crosses), the nucleon (stars) 
and the (filled triangles) 
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FIG. 11: Comparison of the charge density distribution to the 
matter distribution for the unquenched theory at k = 0.156. 
The notation is the same as that of Fig. |7| 



deformation was obtained If the pion cloud makes 
a significant contribution to the charge deformation then 
an analysis with lighter dynamical quarks should increase 
the charge deformation but it should also yield a defor- 
mation in the matter distribution. Clearly, in view of 
such important phenomenological implications, it is im- 
perative to improve the lattice results, in particular for 
the baryons, by using lighter quarks on a larger volume 
both in the quenched and in the unquenched theory, and 
determine the amount of deformation in the chiral limit. 



In order to investigate the importance of dynamical 
quarks, we analyse a subset of the SESAM ^ configura- 
tions. The lattice spacing determined from the rho mass 
in the chiral limit using the SESAM configurations is 
a^^ ~ 2.3 GeV 8] which is the same as for the quenched 
theory at /? = 6.0, and therefore the physical volume is 
the same as in our quenched calculation. As we have 
already explained, we use this determination of the lat- 
tice spacing, which is applicable in both the quenched 
and the unquenched theory, to compare quenched and 
unquenched results. For each of the the two values of the 
hopping parameter, k ~ 0.156 and 0.157, we analysed 100 
SESAM configurations. The ratio of the pion mass to rho 
mass is 0.83 ai k — 0.156 and 0.76 at k = 0.157. These 
values are close to the quenched mass ratios measured 
at K = 0.153 (0.84) and k = 0.154 (0.78) respectively, 
allowing us to make pairwise quenched-unquenched com- 
parisons. Such comparison was made for the charge dis- 
tribution in ref. 4]. In Fig. ^| we compare the mat- 
ter distribution for the quenched [k, — 0.153) and the 
unquenched (k = 0.156) theories. We observe no un- 
quenching effects for the baryons, whereas the matter 
density distribution for the pion and rho increases at the 
origin. Beyond about 0.7 fm the statistical errors, not 
shown in the figure for clarity, become large especially 
for the nucleon and the A"*" . In Fig. ^] we compare the 
unquenched charge and matter density distributions at 
K = 0.156. Similar results with larger statistical errors 
are obtained for the lighter quark mass (k = 0.157), indi- 
cating that the mass dependence of these results is very 
weak. As in the quenched case we observe a faster fall 
off of the matter density distribution as compared to the 
charge density distribution. Whereas unquenching in- 
creases the rho asymmetry for the charge correlator (see 
Fig. 17 of ref. 3), it has no effect on the matter density 
distribution: Fig. [SJj remains unchanged. 



V. PSEUDOSCALAR DENSITY 

The pseudoscalar density is of theoretical interest, es- 
pecially because it can serve as an additional observ- 
able in our comparison with bag model predictions. The 
lattice results in the quenched theory for the pion, the 
rho, the nucleon and the A"*" are shown in Fig. E| for 
m^/mp = 0.84, 0.78 and 0.70. Again we observe only a 
weak dependence on the quark mass. To compare with 
bag model predictions, which will be discussed in the 
following section, it is best to consider the pseudoscalar 
density weighted with . This is shown for both the 
quenched and the unquenched theory at k = 0.153 and 
K = 0.156 respectively in Fig. El A reasonable fit is ob- 
tained using an exponential times a polynomial ansatz. 
We note here that in both the quenched and unquenched 
theory the long tail of the data and the integral of the fit- 
ted ansatz both favor a non-zero integral. In Fig. 1141 we 
show the unquenched pseudoscalar density for the four 
hadrons under consideration normalized to unity at the 
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FIG. 12: The pseudoscalar density ai k — 0.153, 0.154 and 
0.155 for (a) the pion, (b) the rho, (c) the nucleon and the 
(d) A+. 



origin. We observe that the results are very similar for all 
hadrons, with almost no dependence on the quark mass. 
As in the case of the matter distribution dynamical quark 
effects are small for these values of quark masses. Before 
leaving this section we must stress that these results on 
the pseudoscalar density are most useful in showing qual- 
itative features. Since we have used an unimproved Dirac 
operator which has order a chiral symmetry violations, 
we expect sizable ultraviolet corrections. 



VI. COMPARISON OF LATTICE AND BAG 
MODEL RESULTS 

It is interesting to compare our lattice results to those 
obtained in the bag model. We will consider only the 
lowest mode of the free Dirac field in a spherical bag of 
radius R. The radius R is chosen so as to minimize the 
mass, M(R), of the hadron under consideration. There 
are four contributions to the mass: the volume term 
= AttBR^ /3, the zero point energy, Eq = —Zq/R, the 
kinetic energy of the quarks, Eq — riqLu, and the chromo- 
magnetic hyper-fine interaction energy calculated to first 
order in perturbation theory which is proportional to 
o^bag [f^ - is the number of quarks in the hadron, 
ttbag = 5'^/47r where g is the strong coupling constant, 
and LO is the frequency of the lowest mode given ex- 
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FIG. 13: The quenched (left) and unquenched (right) pseu- 
doscalar density at k = 0.153 and k — 0.156 respectively with 
fits to (a -I- bx^ + cx*)e~'"^ (solid line) and bag model results 
(dashed line). 



plicitly below. We use two procedures to fix the bag 
model parameters B, Zq and abag: In the first proce- 
dure, which we will refer to as procedure A, we use the 
lattice values for the masses of the rho, the nucleon and 
the A to fix the three parameters using as an input the 
naive quark mass, niq = l/2a(l/K — 1/kc). In Tabled 
we give the bag parameters determined from this proce- 
dure. In the quenched the quark mass decreases 
the bag model parameters nicely approach the param- 
eters determined using the experimental values of the 
LO meson, the nucleon and the A and mg = also in- 
cluded in the last row of Table HJ For the unquenched 
case, this trend is only clear for the two smallest k val- 
ues (heaviest quarks). For the larger k values this is no 
longer clear and a possible explanation is finite volume 
effects known to become important at these sma ll q uark 
masses. A second procedure suggested in ref. [ig: re- 
ferred to here as procedure B, is to use the bag model 
parameters as obtained in the last row of Table and 
with these parameters fixed, adjust the quark mass so 
that the lattice results for the mass of the nucleon is 
reproduced at the given k value. Since we use the nu- 
cleon mass to define the quark mass in this way, it is 
natural for this discussion to use the nucleon mass to 
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set the lattice scale a. This gives a^^ — 2.04(2) GeV 
for the quenched case and = 1.88(7) for the un- 
quenched theory. The quark mass obtained with proce- 
dure B is denoted by m^^^. Both the bag model and 
the lattice phenomenology are consistent with a linear 
dependence of the nucleon mass on the quark mass in 
the regime we explore. Therefore we expect niq cx m^^^. 
In Fig. El we display the values obtained for m^^s ver- 
sus the naive quark mass. The quenched data nicely fall 
on a straight hne confirming this expectation. We find 
that m^^s is about 2.1 times the naive quark mass. Ex- 
trapolating to niq = we obtain to^'^s _ g 2 MeV 
reasonably close to the expected zero value, especially 
since the error does not include uncertainties of the or- 
der of 10% in setting the lattice scale. The same anal- 
ysis can also be done for the less accurate unquenched 



data. We find m^^s ~ 3to,. Using the SESAM resuhs for 
the nucleon mass at k = 0.156, 0.1565, 0.157 and 0.1575 
and a linear extrapolation to the chiral limit, we obtain 
771^^8 = 31 ± 24 MeV at mq = 0, again close to zero. In 
the same figure we also show the bag model results for the 
matter and pseudoscalar densities calculated with proce- 
dure A using the parameters extracted from fitting the 
quenched data at k = 0.153, 0.154 and k = 0.155. The 
bag model results show a stronger dependence on the 
quark mass in comparison with the corresponding lattice 
results discussed in the previous sections. 

In order to evaluate the various density distributions 
we expand the quark fields in terms of bag eigenmodes 
keeping only the lowest mode [Tsf : 



q{r) = ^(6^V'.(r) + rf:^</'.(r) 



TABLE 1: Bag model parameters extracted from fitting the 
masses of the rho, the nucleon and the A. The bag radii, Rtt, 
Rp, Rn and Ra for the pion, the rho, the nucleon and the A 
are given in GeV~^ respectively. 



Quenched: /3 = 6.0, 16^ x 32 


(GeV) Bi/^(GeV) Zo abag 


Rtv Rp Rn Ra 


0.174 0.209 2.274 0.169 
0.131 0.182 2.245 0.196 
0.088 0.169 2.228 0.239 


2.71 2.99 3.48 3.59 
3.09 3.47 4.00 4.15 
3.26 3.76 4.29 4.49 


Nf =2: 13 = 5.6, 16^ x 32 


m, (GeV) Bi/''(GeV) Zq abag 


Rn Rp Rn Ra 


0.096 0.210 2.120 0.110 
0.076 0.195 2.096 0.136 
0.057 0.201 2.383 0.188 
0.038 0.171 2.227 0.188 


2.86 3.02 3.53 3.60 
3.06 3.28 3.80 3.89 
2.74 3.09 3.60 3.73 
3.33 3.70 4.28 4.43 


Continuum 


m, (GeV) Bi/4(GeV) Zo abag 


Rn Rp Rn Ra 


0.146 1.86 0.55 


3.33 4.67 4.97 5.44 



fir) = 



9{r) = 



« I{r)Us 
-g{r)(j.rUs 

N{x) 



0,(r) = C7o^:(r) 



'Air 
N{x) 



1/2 



jo{xr/R) 



1/2 



ji (xr/ R) 



(7) 



where Us are two component Pauli spinors, C = 7270 is 
the charge operator, ji are spherical Bessel functions, m 
is the value of the quark mass which here is cither or 
TO^^s a,nd N{x) is fixed by normalizing the cigcnmodc in 
the bag. The frequency, uj{m,R), of the lowest mode is 



given hy uj — l/R {x 



2^2-ji/2 .yjfj^gj.g x{mR) is the 



solution of the eigenvalue equation (X — mR—ojR) tan a; = 
X. The superscript c on the quark annihilation operators 
hi and on the anti-quark creation operators ^ denotes 
flavour and colour quantum numbers. In terms of the 
lowest eigenmode the density operators are given by 



p,„(r) = 5](6^t5c 



/575(r) 



d^d^) (fir)' 



9{rf) 



(8) 
(9) 



-dlUl.Ula^.vUs') f{r)g{r) (10) 

and the charge, matter and pseudoscalar correlators by 
the following expressions: 
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3lM)3iM') \P 
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{fir 



9{ry 



V 1 / 



{f{r'f+g{r'?) 

(11) 
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FIG. 15: Upper: The quark mass, mj^^^, needed to obtain the 
lattice result for the nucleon mass versus the naive quark mass 
ruq. The best straight line fits to both the quenched (dashed) 
and unquenched (dotted) case are also displayed. Lower: The 
mass dependence of the matter (left) and pseudoscalar (right) 
density distributions for the nucleon within the bag model 
with parameters fixed using procedure A at k = 0.153, 0.154 
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4f.f7(r)3(r) 
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(13) 

with constants C, C and C" independent of the hadron 
state. The minus sign for mesons in the charge correlator 
arises because in the density operator the term involving 
anti-quarks comes with a negative sign. This term is pos- 
itive in the case of the matter density giving the same sign 
for mesons and baryons. This shows exphcitly that the 
difference between charge and matter density distribu- 



tions is due to the opposite sign of the lower components 
of the Dirac spinor which it is a relativistic effect. 

In the case of the pseudoscalar correlator given by 
Eq. 1131 the results were obtained by averaging over spin 
projections of the physical states i.e. for the rho wc av- 
eraged over the states with = ±1 and ~ Q- In 
order to understand the signs for the pseudoscalar den- 
sity consider the pion where the spin of the quark and the 
anti-quark are opposite. Since |7r) ~ l/\/2^^(6" ^dl^ — 
6" ^(i'J_^)|0), we have from Eq. ^| a contribution of the 
form ±a.TU±U:fa* .t'U- yielding the result — f.f' and 
thus an overall positive sign. The rho state with = 
is orthogonal to the pion resulting in the opposite sign. 
Similar considerations lead to the opposite sign between 
the nucleon and the A"*" . 

In Fig. El we compare the lattice results for the charge 
and matter density distributions with the bag model re- 
sults from Eqs. ^2 and ^1 The dashed lines show the 
results obtained with procedure A. The dotted lines are 
the bag model results obtained within procedure B. Both 
procedures fail to reproduce the correct radial depen- 
dence of the charge distribution. The results obtained 
using procedure A provide an overall better description 
in the case of the baryon matter distribution. To make 
the comparison quantitative, we evaluate the root mean 
square radius, which provides a measure of the width of 
the distributions. In the bag model one usually computes 
the expectation value of for the lowest eigenmode. In 
our previous study of the charge density distribution 
we used the quark model definition for the charge radius 
which for mesons is given by 



^eg((r, -Rem)') 
E.e.Jd'r (r/2)^C,„(r) 



(14) 



where Rem is the coordinate of the center of mass and 
Bq is the electric charge of the quarks. A correspond- 
ing definition for baryons can only be used if one knows 
the charge density distribution in terms of the two rel- 
ative coordinates which requires the evaluation of three 
current correlation functions. Here we only evaluate the 
one-density baryon charge distribution and therefore the 
quark model definition can not be applied. For simplicity 
and for direct comparison with the the bag model radial 
width we calculate for both mesons and baryons 



r- = - 



iJjfrvHJrir) 
2 Jd^rCrir) 



(15) 



where for F = 70 we obtain the charge mean square ra- 
dius, r^^, and for F = / the matter mean square radius, 
r^„. For degenerate quarks and for the case where the 
meson wave function is a product of single particle radial 
wave functions like in the bag- model, the factor of 1/2 
corrects for summing over the charge root mean square 
radius of each quark as it is done in Eq. 1151 in using the 
relative quark distance square. 
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The lattice results for the charge and matter radii 
are collected in Tabled The main observation is that 
hadron sizes show very Httle quark mass dependence in 
the quenched case, at least for the range of quark masses 
considered, but increase under the effect of dynamical 
quarks. One should keep in mind, however, that finite 
volume effects on these values have not been investigated 
here, but could be significant at the smaller quark masses. 

Despite the failure of the bag model in describing the 
individual radial shape of the distributions it produces 
reasonable results for the relative widths of the charge 
to the matter distribution as can be seen from Table ITTll 
Both lattice and bag model consistently predict a broader 
charge than matter distribution, with a very weak mass 
dependence. In Table lllll we also give the bag-model 
results for the charge radius which show a stronger de- 
pendence on the quark mass as well as on the channel 
than the lattice results. 
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FIG. 16: Comparison of quenched charge (left) and matter 
distribution (right) at k = 0.153 with bag model results. The 
dashed line is obtained by using procedure A as described in 
the text and the dotted line by using procedure B. The solid 
line is a fit to aexp(— mr"). 

In Fig. El in addition to the bag-model results, we 
have also included a fit of the lattice data to the ansatz 
exp(— mr"). As it can be seen this simple ansatz pro- 
vides the best description to the data. In fact, except 
for the pion, a good fit to all other charge and matter 



correlators is obtained by taking a = 1, which gives a 
per degree of freedom < 1. In this case we can iden- 
tify asymptotically the mass m of the pure exponential 
to the mass of a propagating meson. To fit the matter 
distributions we find that we need a mass of approxi- 
mately twice that required to fit the charge distributions 
of these hadrons. The asymptotic behaviour of the den- 
sity correlators can be analysed using a tree graph ap- 
proximation based on the fact that at large distances 
light hadrons dominate . In our work the ratio of the 
pion to the rho mass is alwa ys l arger than one-half. In 
this case, the analysis of ref. [13 shows that the density 
correlators should decay exponentially with the mass of 
the rho meson, expect for the rho correlator whose expo- 
nential decay at very large distances is governed by the 
pion mass. A similar analysis for the matter correlator 
predicts asymptotically an exponential decay determined 
from the lightest scalar in the theory. In the continuum 
the lightest scalar (/o) is heavier than the rho by 200 MeV 
and one thus expects asymptotically a faster fall off of the 
matter correlator as compared to the charge correlator. 
Unfortunately, the values we obtain for the masses, to, 
of the exponential fall off from the fits can not be iden- 
tified with the mesons of the theory. In general we find 
slower exponential decay than expected from the mass of 
the lightest meson which should dominate the asymptotic 
behaviour. Presumably, our lattice is not large enough to 
probe the asymptotic behaviour of the correlators. The 
pion correlators require a < 1, or the sum of two ex- 
ponentials governing the short and long distance decays 
respectively. If we nevertheless use the single exponen- 
tial ansatz with a = 1 also for this channel the ratio of 
mass values needed to fit the charge and matter density 
distributions is about 1.4 instead of two as for the other 
hadrons decreasing slightly with the quark mass. This 
value of ~ 1.4 is thus in accord with the value of ~ 1.6 
found from the study of heavy light mesons in ref. [T^ . 

In the case of the pseudoscalar density the bag model 
predicts that the integral J d?r C^(r) is zero. This pre- 
diction of the bag model is clearly seen in Fig.^] It was 
for this reason that in Fig.^jwe chose to show the lattice 
quenched and unquenched results for the pseudoscalar 
density weighted by r^. As can be seen from Fig. ^1 for 
both the quenched and the unquenched results, the bag 
model prediction is worse for the pion, since the lattice 
correlator turns negative at a larger distance. This is not 
surprising since the bag model is known to be worse for 
the pion. For the baryons on the other hand the node 
of the bag model distribution coincides with that of the 
lattice data. In all cases the lattice data have a long tail 
which is not reproduced in the bag model, and which 
favors a non-zero integral. However a careful thermody- 
namic and continuum extrapolation, on larger and finer 
lattices, is required for this quantity, especially with Wil- 
son fermions as used here. For the smaller quark masses 
noise at the tail of the distribution makes the situation 
even worse. 



12 



TABLE II: Square root of the charge (r^h) a-nd matter ( 
r^) radii in fm extracted from the lattice data. We use the 
nucleon mass to convert to physical units. 



Quenched: /? = 6.0, 16'^ x 32 



K 


0.153 

^FTh ^ 


0.154 

\/^ 


0.155 


TV 
P 

n 

A+ 


0.444(4) 0.396(6) 
0.457(5) 0.397(18) 
0.457(4) 0.383(16) 
0.464(6) 0.390(19) 


0.454(4) 0.417(6) 
0.470(6) 0.404(37) 
0.465(6) 0.384(31) 
0.469(9) 0.360(76) 


0.465(8) 0.440(7) 
0.482(11) 
0.473(11) 
0.447(77) 


Nf = 2: 13 = 5.6, 16^ x 32 


K 


0.156 


0.157 

^/rlh 




TV 
P 

n 

A+ 


0.471(3) 0.411(7) 
0.481(6) 0.370(52) 
0.485(5) 0.427(14) 
0.489(9) 0.384(55) 


0.490(3) 0.439(7) 
0.510(6) 0.373(121) 
0.509(7) 0.414(57) 
0.522(9) 0.436(102) 





TABLE III: Ratio of the expectation value of for the charge 
distribution over that for the matter distribution in quenched 
lattice QCD and bag model using procedure A. For «: = 0.155 
we give the ratio only for the pion since for the other hadrons 
the matter radius was too noisy. The root mean square charge 
radius in fm in the bag-model using procedure A is also given. 





K 


0.153 




0.154 


0.155 






lattice 


bag 


lattice bag 


lattice 


bag 


TV 


1.24(3) 


1.35 


1.19(2) 1.36 


1.10(2) 


1.37 


P 


1.33(7) 


1.35 


1.35(17) 1.36 




1.37 


n 


1.40(8) 


1.34 


1.45(16) 1.35 




1.36 


A 


1.41(9) 


1.34 


1.7(5) 1.35 




1.36 






/ ' ch 


in the bag model (fm) 




K 


0.153 


0.154 


0.155 


TV 


0.383 




0.439 


0.466 




P 


0.422 




0.492 


0.537 




n 


0.488 




0.565 


0.612 




A 


0.503 




0.586 


0.639 





VII. CONCLUSIONS 

To check our lattice procedure we performed an analy- 
sis of the charge density distributions using various types 
of sources. For a lattice of temporal extent 32a consid- 
ered here and within our statistics we found that local 
and Wuppertal smeared sources give consistent results 
when we insert the current operator at time separation 
from the source of f = 8a, which is the maximal allowed 



separation given the anti-periodic boundary conditions in 
the temporal direction. The results are also in agreement 
with those obtained using Dirichlet boundary conditions 
which allow larger time separations. Therefore, for the 
lattice parameters used in this study, local sources are a 
suitable choice for the evaluation of density distributions 
since they have less gauge noise than smeared ones and 
the temporal extent is large enough to make contribu- 
tions from high excited states negligible. Wall sources 
are shown to filter the ground state less effectively and 
at the maximal time separation of 8a yield results that 
are not in agreement with the rest. Using the lattice 
conserved current, which has a different a-dependence 
as compared to the continuum current, these correlators 
show finite lattice spacing effects near the origin. 

Our main conclusion from the comparison of quenched 
and unquenched results for the charge, the matter and 
the pseudoscalar densities is that no sizable unquench- 
ing effects are observed for quark masses in the range of 
300-100 MeV. The charge density distribution is, in all 
cases, broader than the matter density. This is what is 
expected if a tree level classical approximation is used to 
describe the asymptotic behaviour of these correlators. 
However the masses of the propagating lightest mesons 
which determine the exponential fall off can not be iden- 
tified with the mesons of the theory. To observe the true 
asymptotic behaviour larger lattices are needed. For the 
pion, the lattice results are in agreement with the ex- 
perimentally extracted ratio of the charge to the matter 
radius of 1.15(9) For baryons, the lattice indicates a 
charge radius about 20% larger than the matter radius. 
This effect is well reproduced by the bag model. However 
the bag model does a poor job in the description of the 
radial dependence, especially in the case of the mesons. 
Instead, the charge and the matter distributions are well 
described by the simple ansatz exp(— mr"), with a — 1 
except for the pion. One prediction of the bag model is 
that the the volume integral of the pseudoscalar density 
is zero. Lattice data for both the quenched and the un- 
quenched theory favor a negative value. A larger lattice 
is required to settle this issue. 

The deformation seen in the rho charge distribution is 
absent in the matter distribution, both in the quenched 
and the unquenched theory. This observation suggests 
a relativistic origin for the deformation. This important 
issue deserves a more extended study, with lighter quarks 
and larger volumes. 
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